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We present a mechanism for a generic and powerful force of as¬ 
sembly and mobility for transmembrane proteins in lipid bilayers. 
This force is a pre-transition (or pre-melting) effect for the first- 
order transition between ordered and disordered phases in the host 
membrane. Using large scale molecular simulation, we show that a 
protein with hydrophobic thickness equal to that of the disordered 
phase embedded in an ordered bilayer stabilizes a microscopic or¬ 
der-disorder interface, and the stiffness of that interface is finite. 
When two such proteins approach each other, they assemble be¬ 
cause assembly reduces the net interfacial free energy. In analogy 
with the hydrophobic effect, we refer to this phenomenon as the 
“orderphobic effect.” The effect is mediated by proximity to the 
order-disorder phase transition and the size and hydrophobic mis¬ 
match of the protein. The strength and range of forces arising from 
the orderphobic effect are significantly larger than those that could 
arise from membrane elasticity for the membranes we examine. 

lipid bilayers, phase transition, hydrophobic mismatch, orderphobe 

Abbreviations: DPPC: dipalmitoyl phosphatidylcholine 


This paper presents implications of the first-order or¬ 
der-disorder phase transition in lipid bilayers. The fluid mo¬ 
saic model [1] and the lipid raft hypothesis [2] have guided 
intuition on how proteins diffuse in biological membranes — 
ordered clusters floating in an otherwise disordered fluid mem¬ 
brane [3,4]. However, recent advances show that the state 
of membranes containing transmembrane proteins is ordered, 
even gel-like [5-10]. How then do transmembrane proteins dif¬ 
fuse and assemble within this relatively rigid material? Here, 
we argue that this question is answered by the fact that a 
transmembrane protein in an ordered bilayer can induce ef¬ 
fects that resemble pre-melting [11-13]. Specifically, with 
molecular simulation, we show that within an otherwise or¬ 
dered membrane phase, mesoscopic disordered domains sur¬ 
round proteins that favor disordered states. We find, impor¬ 
tantly, that the boundary of the domains resembles a sta¬ 
ble, fluctuating order-disorder interface. The dynamic equi¬ 
librium established at the boundary allows the protein and its 
surrounding domain to diffuse. 

Moreover, because the interface has a finite stiffness, 
neighboring proteins can experience a membrane-induced 
force of adhesion, an attractive force that is distinctly stronger 
and can act over significantly larger lengths than those that 
can arise from simple elastic deformations of the membrane 
[14-18]. This force between transmembrane proteins is anal¬ 
ogous to forces of interaction between hydrated hydrophobic 
objects. In particular, extended hydrophobic surfaces in wa¬ 
ter can nucleate vapor-liquid-like interfaces. In the presence 
of such interfaces, hydrophobic objects cluster to reduce the 
net interfacial free energy. This microscopic pre-transition ef¬ 
fect manifesting the liquid-vapor phase transition can occur 
at ambient conditions [19-26]. 

In the transmembrane case, we show here that a protein 
favoring the disordered phase creates a similar pre-transition 
effect. In this case it manifests the order-disorder transition 


of a lipid bilayer. Like the raft hypothesis, therefore, clusters 
do indeed form, but the mechanism for their assembly and 
mobility emerge as consequences of order-disorder interfaces 
in an otherwise ordered phase. We refer to this phenomenon 
as the “orderphobic effect.” 


The order-disorder transition is a first-order phase transition. 

We choose the MARTINI model of hydrated dipalmitoyl phos¬ 
phatidylcholine (DPPC) lipid bilayers [27] to illustrate the or¬ 
derphobic effect. See Methods. This membrane model exhibits 
an ordered phase and a disordered phase. Fig. 1A contrasts 
configurations from the two phases, and it shows our esti¬ 
mated phase boundary between the two phases. The ordered 
phase has regular tail packing compared to the disorganized 
tail arrangement of the disordered phase. A consequence of 
the regular tail packing is that hydrophobic thickness of the 
ordered phase, V Q , is larger than that of the disordered phase, 
V a. Correspondingly, the area per lipid in the ordered phase 
is smaller than that in the disordered phase. 

Rendering the end particles of all the lipid chains in one 
of the two monolayers provides a convenient visual represen¬ 
tation that distinguishes the two phases. These tail-end par¬ 
ticles appear hexagonally-packed in the ordered phase and 
randomly arranged in the disordered phase. Regions that ap¬ 
pear empty in this rendering are in fact typically filled by non 
tail-end particles or by tail-end particles from the other lipid 
monolayer. 


Significance 

The clustering of proteins in biological membranes is a con¬ 
trolling factor in processes such as endo- and exo-cytosis, cell 
signaling, and immunological synapses. Yet physical principles 
governing this organization are incompletely understood. Here, 
we address some of this uncertainty by demonstrating a gen¬ 
eral mechanism for mobility and powerful forces of assembly for 
transmembrane proteins. 
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Fig. 1 . (A) Order-disorder phase diagram in the tension-temperature, A-T, plane. 
The lateral pressure across the membrane is —A. Points are estimated from 10 inde¬ 
pendent heating runs like those illustrated in Fig. SI for a periodic system with 128 
lipids. Insets are cross sections showing configurations of a bilayer with 3200 lipids 
in the ordered and disordered phases. The heads are colored gray while the tails are 
colored pink. Water particles are omitted for clarity. The hydrophobic thicknesses, 
T> 0 and T>d’ are the average vertical distances from the first tail particle of the upper 
monolayer to that of the lower monolayer. A macroscopic membrane buckles for all 
A < 0. Snapshots of the last tail beads in one monolayer of each phase are shown 
to illustrate the difference in packing. (B) Snapshot of a system showing coexistence 
between the ordered and disordered phases. The gray contour line indicates the loca¬ 
tion of the interface separating the ordered and disordered regions. The snapshot is 
a top view of the bilayer showing the tail-end particles of each lipid in one monolayer. 
(C) Fourier spectrum of the fluctuations of the instantaneous order-disorder interface. 
The line is the small-fc capillarity-theory behavior with 7 = 11.5 pN. 

To quantify the distinctions between the two phases, 
we consider a local rotational-invariant [28-30], fa = 

I ( 1 / 6 ) Syenn(z) exp | 2 , where Oij is the angle between 
an arbitrary axis and a vector connecting tail-end particle l to 
tail-end particle j, and the summation is over the six nearest 
neighbors of particle l. The equilibrium average, (fa), is 1 for 
a perfect hexagonal packing, and it is 1/6 or smaller in the 
absence of bond-orientation correlations. Small periodically 
replicated samples of the DPPC hydrated membrane exhibit 
hysteretic changes in area per lipid and in (fa) during heating 
and cooling. See Supporting Information (SI), and Refs. [31] 


and [32]. To establish whether the first-order-like behavior 
persists to large scales and thus actually manifests a phase 
transition, we consider larger systems and the behavior of the 
interface that separates the ordered and disordered phases. 

Fig. IB shows coexistence for a system size of IV = 3900 
lipids with an interface between the two phases. To analyze 
interfacial fluctuations, we first identify the location of the 
interface at each instant. This location is found with a two- 
dimensional version of the three-dimensional constructions de¬ 
scribed in Refs. [13] and [33]. Specifically, and as discussed in 
Methods , the interface is the line in the plane of the bilayer 
with an intermediate coarse-grained value of the orientational- 
order density, 

<M r ) = 53<M(r- r i)- [1] 

l 

where vi is the position of the Ith tail-end particle projected 
onto a plane parallel to that of the bilayer, r is a two- 
dimensional vector specifying a position in that plane, and 
£(r) is Dirac’s delta function. We focus on this field rather 
than the tail-end number density, p(r) = ^ 5(r — rj), be¬ 
cause the difference between the two phases is larger for typ¬ 
ical orientational-order than for typical tail-end density. 

Fig. 1C shows the Fourier spectrum of the height fluctua¬ 
tions of this interface, (\6hk\ 2 ). Two different system sizes are 
studied, with the larger system having approximately dou¬ 
ble the interface length of the smaller system. The Fourier 
component 5hk is related to the height fluctuation 5h(x) as 
5h(x) = 5hk exp (ikx) where x is a point along the hor¬ 
izontal in Fig. IB. Here, 0 ^ x ^ L, and L is the box 
length. With periodic boundary conditions, k = 2r:m/L , 
m = 0, d=l, ± 2 , • • •. According to capillarity theory for crys¬ 
tal-liquid interfaces [34,35], (\Shk\ 2 ) ~ feT/Ly/c 2 for small k, 
with &b being the Boltzmann’s constant. 

Given the proportionality with l/k 2 at small k (i.e., wave¬ 
lengths larger than 10 nm), comparison of the proportionality 
constants from simulation and capillarity theory determines 
the interfacial stiffness, yielding 7 = 11.5 zb 0.46 pN. This 
value is significantly larger than the prior estimate of interfa¬ 
cial stiffness for this model, 3d=2pN [31]. That prior estimate 
was obtained from simulations of coarsening of the ordered 
phase. 

Because the ordered phase has a hexagonal packing, the 
interfacial stiffness depends on the angle between the inter¬ 
face and the lattice of the ordered phase. For a hexagonal 
lattice, there are three symmetric orientations for which the 
interfacial stiffnesses are equal. We will see that for the model 
we have simulated there appears to be only little angle depen¬ 
dence. Irrespective of that angle dependence, the stability of 
the interface and the quantitative consistency with capillary 
scaling provide our evidence for the order-disorder transition 
being a first-order transition in the model we have simulated. 

The system sizes we have considered contain up to 10 7 
particles, allowing for membranes with N ~ 10 4 lipids, and re¬ 
quiring 10 /is to equilibrate. As such, our straightforward sim¬ 
ulations are unable to determine whether the ordered phase 
is hexatic or crystal because correlation functions that would 
distinguish one from the other [36] require equilibrating sys¬ 
tems at least 10 times larger [37]. Similarly, we are unable to 
determine the range of conditions for which the membranes or¬ 
ganize with ripples and with tilted lipids [38,39]. Presumably, 
the ordered domain of the phase diagram in Fig. 1A, parti¬ 
tions into several subdomains coinciding with one or more of 
these possibilities. With advanced sampling techniques [40], 
free energy functions of characteristic order parameters can 
be computed to estimate the positions of boundaries between 
these various ordered behaviors. Here, we do not pursue this 
additional level of detail in the phase diagram because the ad- 
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ditional boundaries refer to continuous transitions [38]. It is 
only the first-order transition, with its discontinuous change 
between ordered and disordered phases, that supports coexis¬ 
tence with a finite interfacial stiffness, and it is this stiffness 
that results in the orderphobic effect, which we turn to now. 



Protein 




r—R [nm] 

Fig. 2. Model proteins in the bilayer. (A) Idealized cylindrical protein-like so¬ 
lutes with radius R and hydrophobic thickness I (magenta). The hydrophilic caps 
of the protein are shown in white. (B) Cross section of the lipid bilayer in the or¬ 
dered phase containing a model protein of radius 2.7 nm with a hydrophobic thickness 
.£=2.3 nm < T>d- (C) The radial variation of the order parameters ((f) (r)) (right 
axis) and (</>(r))/(p(r)) (left axis) show disorder in the vicinity of the protein of 
radius 1.9 nm. (D) Comparison of the radial order parameter variation for three dif¬ 
ferent proteins shows an increase in the extent of the induced disorder region with 
protein radius. 


Transmembrane proteins can disfavor the ordered membrane. 

A disordering (i.e., orderphobic) transmembrane protein is 
one that solvates more favorably in the disordered phase than 
in the ordered phase. The disordering effect of the protein 
could be produced by specific side chain structures. See SI. 
Here, in the main text, we consider a simpler mechanism. In 
particular, we have chosen to focus on the size of the protein’s 
hydrophobic thickness and the extent to which that thick¬ 
ness matches the thickness of the membrane’s hydrophobic 
layer [41,42]. See Fig. 2. 

The membrane’s hydrophobic layer is thicker in the or¬ 
dered state than in the disordered state. For instance, at zero 
lateral pressure and 294 K in the model DPPC membrane, 
we find that the average thicknesses of the hydrophobic lay¬ 
ers in the ordered and disordered states are V Q m 3.1 nm and 
Vd — 2.6 nm, respectively. A transmembrane protein with 
hydrophobic thickness of size I ~ 2.6 nm will therefore favor 
the structure of the disordered phase. If the protein is large 


enough, it can melt the ordered phase near the protein and 
result in the formation of an order-disorder interface. 

Spatial variation of the order parameter field character¬ 
izes the spatial extent of the pre-melting layer. To evaluate 
whether a model protein is nucleating a disordered domain 
in its vicinity, we calculate the average of the orientational- 
order density held as a function of r — |r|, (0(r)) (right axis 
of Fig. 2C). It exhibits oscillations manifesting the atom¬ 
istic granularity of the system. Dividing by the mean density, 
(p(r)), largely removes these oscillations. 

A profile of this ratio in the vicinity of the protein is de¬ 
picted in Fig. 2C (left axis). It changes approximately sig¬ 
moidally, connecting its values of 0.15 and 0.45 in the disor¬ 
dered and ordered phases, respectively. The shape of the pro¬ 
file suggests the formation of an order-disorder interface [43]. 
Further, the increase in the spatial extent of the disordered 
region with the increasing size of the protein Fig. 2D is in¬ 
dicative of length scale dependent broadening effects brought 
about by capillary fluctuations. These impressions can be 
quantified by analyzing fluctuations of the instantaneous in¬ 
terface, which we turn to now. 

An orderphobic protein nucleates a fluctuating order-disorder 
interface. Fig. 3A shows a configuration of the instantaneous 
interface that forms around the orderphobic protein shown in 
Fig. 2B. The interface is identified as described above. A 
video of its dynamics is provided at https ://goo .gl/NBQJP9. 




Fig. 3. Soft order-disorder interface: (A) Arrangement of the tail end particles 
of the top monolayer corresponding to the protein in Fig: 2B. Far away from the 
protein, the tail end particles show hexagonal-like packing and are in the ordered 
state. Proximal to the protein, it can be seen that the tail end particles are ran¬ 
domly arranged, and resemble the disordered phase. The line connected by the black 
points denotes the instantaneous order-disorder interface. (B) The fluctuations in 
the radius of the order-disorder interface are consistent with the fluctuations of a free 
order-disorder interface at coexistence. IZo is the mean radius of the order-disorder 
interface surrounding a model protein of radius R. 
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Fig. 4. Demonstration of the orderphobic force: Two proteins separated by a center-to-center distance of 14 nm are simulated at 309 K. Snapshots at various times reveal 
the process of assembly in which the two order-disorder interfaces merge into a single interface. 


As is common in crystal-liquid interfaces, the interface nu¬ 
cleated by an orderphobic protein may exhibit hexagonal 
faceting [34], remnants of which can be observed in Fig. 3A. 

The mean interface is a circle of radius IZq. Fourier 
analysis of fluctuations about that circle yields a spectrum 
of components. To the extent that these fluctuations obey 
statistics of capillary wave theory for circular interface, the 
mean-square fluctuation for the kth component is {\5IZk \ 2 ) — 
k-BT/2ir^k 2 'JZo, where k — m/IZo and m m ±1, ± 2 , • • •, and 
7 is the order-disorder interfacial stiffness, neglecting the de¬ 
pendence on the angle between the interface and the lattice. 
The discrete values of k reflect periodic boundary conditions 
going full circle around the model protein. 

In Fig. 3B, we use the interfacial stiffness from the free 
interface (7 = 11.5 pN) separating coexisting ordered and dis¬ 
ordered phases with the capillary theory expression, and its 
corresponding spectrum, to compare with the spectrum of the 
protein-induced interface. The agreement between the theory, 
free interface and the protein-induced interface is good, and 
it improves as the radius of the orderphobic protein increases 
and the wave vector k decreases. This agreement indicates 
that the orderphobic protein does indeed nucleate an interface 
manifesting the order-disorder transition. The deviations of 
the fluctuations of the free interface from capillary wave the¬ 
ory occur for k > 0.8 nm -1 , corresponding to wavelengths 
2?r//c < 7 nm, and a mean interface radius IZo < 1 nm. In¬ 
deed, Fig. 2 suggests that even a small protein of radius 0.5 
nm, which supports an interface of radius IZo ~ 1.2 nm, is suf¬ 
ficient to induce an order-disorder interface with fluctuations 
consistent with capillary theory. 

The orderphobic effect generates forces of assembly and fa¬ 
cilitates protein mobility. Fig. 4 shows three snapshots from a 
typical trajectory initiated with two orderphobic proteins of 
radius 1.5 nm separated by a distance of 14 nm. Each induces 
a disordered region in its vicinity, with soft interfaces sepa¬ 
rating the ordered and disordered regions. The free energy of 
the separated state is approximately j(Pi + P 2 ), where Pi is 
the perimeter of the order-disorder interface around protein i. 
On average, {Pi) = 2 tt1Zo- After a few hundred nanoseconds, 
a fluctuation occurs where the two interfaces combine. While 
the single large interface remains intact, the finite tension of 
the interface pulls the two proteins together. Eventually, the 
tension pulls the two proteins together with a final perimeter, 
Pf, that is typically much smaller than Pi + P 2 . A video of 
its dynamics is provided at https://goo.gl/HXS0j7. 

After the separated interfaces join, the assembly process 
occurs on the time scale of microseconds. This time is required 
for the proteins to push away lipids that lie in the path of the 
assembling proteins. Given this time scale, a reversible work 
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calculation of the binding free energy would best control both 
the distance and the number of lipids between the proteins. 
Moreover, the evident role of interfacial fluctuations indicates 
that the transition state ensemble for assembly must involve 
and interplay between inter-protein separations and lipid or¬ 
dering as well as lipid concentration. 

While we leave the study of reversible work surfaces and 
transition state ensembles to future work, it seems already 
clear that the net driving force for assembly is large compared 
to thermal energies. For example, with a model orderphobic 
protein radius of 1.5nm, we find 7 ((Pf) — 2 (Pl)) « — 30/cbT. 
The range over which the force acts is given by the average 
radius of the two interfaces, 21Zo- This range is further ampli¬ 
fied by the width of the interface, which is of G{IZ 0 ) for one¬ 
dimensional interfaces in two-dimensional systems [44]. The 
typical range is ~ 10 to 30 nm. In comparison, given the elas¬ 
tic moduli of the membranes we consider, elastic responses 
will generate attractive forces between transmembrane pro¬ 
teins that are much smaller in strength and range, typically 
—5 /cbT and 1 nm, respectively [18,45]. Moreover, similarly 
weak and short ranged forces are found from solvation theory 
that accounts for linear response in microscopic detail while 
not accounting for the possibility of an underlying phase tran¬ 
sition [46]. 

As in the hydrophobic effect [19], the strength and range 
of the orderphobic force leverages the power of a phase transi¬ 
tion, depending in this case on the ability of the orderphobic 
protein to induce a disordered layer in its vicinity. This ability 
depends upon the proximity to the membrane’s phase tran¬ 
sition, and, for the simple protein models considered in this 
paper, it depends upon the protein’s radius and hydrophobic 
mismatch with the membrane. The spatial extent of the dis¬ 
ordered region increases with proximity to phase coexistence 
as shown in Fig. 5A. 

Furthermore, Fig. 5B shows that the strength of the effect 
is maximal for a hydrophobic thickness equal to that of the 
disordered phase, and it decreases as the hydrophobic thick¬ 
ness approaches that of the ordered phase. In the case of zero 
mismatch, i.e., I — P Q , the value of the order parameter in 
the vicinity of the protein is consistent with that of a pure 
bilayer in the ordered state. Therefore, the model proteins 
with zero mismatch, do not induce a disordered region, and 
the orderphobic effect vanishes. See Figs. 5B and 5D. 

Fig. 4 also shows that the orderphobic effect produces ex¬ 
cess mobility, by proteins melting order in a surrounding mi¬ 
croscopic layer and by facilitating the motions of neighboring 
proteins. This finding explains how protein mobility and reor¬ 
ganization can be relatively facile in the so-called “gel” phases 
of membranes. Further information on this phenomenon is 
provided in SI. Our prediction of enhanced lipid mobility sur- 
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i = 2.3 nm, T = 294 K 


t — 3.1 nm, T = 294 K 


Fig. 5. Strength of the orderphobic force: (A) Radial variation of the order parameter showing the extent of the disordered region as a function of temperature, for a 
protein of radius 1.9 nm and hydrophobic thickness 2.3 nm. The extent of the disordered region increases as the melting temperature is approached, at zero surface tension. 
(B) Comparison of the radial variation of the order parameter for different hydrophobic mismatches. Proteins with no mismatch do not create any disordered region. (C) 
Arrangement of lipids around a protein with negative mismatch. (D) Arrangement of lipids around a protein with zero mismatch. 


rounding orderphobic proteins may be amenable to experi¬ 
mental tests by single molecule tracking techniques [47]. 

Implications of the orderphobic effect and related phenomena 
in biological membranes. Biological membranes and trans¬ 
membrane proteins are far more complicated than the models 
considered in this paper. Part of the complexity is associated 
with multiple components, which can be studied outside of 
biological contexts. For example, we anticipate that the or¬ 
derphobic effect will be useful in understanding the phase be¬ 
havior that results from mixing cholesterol with pure or mul¬ 
ticomponent lipid bilayers [48-51]. In this case, cholesterol 
with small hydrophilic heads and short hydrophobic tails has 
the propensity to induce disorder in the ordered phase [50]. 

Transmembrane proteins in biological membranes have 
structures more complicated than those of our simple model 
proteins. The side chains of these proteins can affect the 
packing of lipid chains. To the extent that lipid packing is 
dirupted, even small ct-helix proteins can be orderphobic. Ev¬ 
idence for this assessment is provided in SI. How associated 
orderphobic effects lead to clustering and related phenomena 
for transmembrane a-helices merit future study. 

Further, there is a dual to the orderphobic effect: A trans¬ 
membrane protein in the disordered phase that favors the or¬ 
dered phase can nucleate an ordered region and order-disorder 
interface. For example, one of our model proteins with a pos¬ 
itive mismatch [I — T> 0 ) would induce order in its vicinity. 
Interfaces separating the ordered and disordered regions will 
again provide a force for assembly. This case corresponds to 
the situation of lipid rafts [2], which consists of ordered do¬ 
mains floating in otherwise disordered membranes. 

Finally, we speculate that the orderphobic effect plays im¬ 
portant roles in membrane fusion and cell signaling [52-55]. In 
the case of fusion, it would appear that one important role is to 
promote fluctuations in an otherwise stable membrane. Oth¬ 
erwise, it is difficult to conceive of a mechanism by which ther¬ 
mal agitation would be sufficient to destabilize microscopic 
sections of membranes. Such destabilization seems necessary 
for initiating and facilitating membrane fusion. Many pro¬ 
teins are involved in such processes [52,56,57], but it may not 
be a coincidence that the hydrophobic thicknesses of SNARE 


proteins are 25% smaller than that of the ordered membrane 
states [58,59]. 

Materials and Methods 

Molecular simulations. We simulate the MARTINI coarse-grained force field using 
the GROM ACS molecular dynamics package [27,60]. 'Antifreeze' particles are added 
to the solvent to ensure that the solvent does not freeze over the temperature range 
considered in the simulations as in Ref. [27]. Thermostats and barostats control tem¬ 
perature and pressure, and checks were performed to assure that different thermostats 
and barostats yielded similar results [40]. The hydrophobic cores of our idealized pro¬ 
teins are constructed using the same coarse-grained beads as the lipid tails (particle 
Cl in the MARTINI topology [27]). Similarly, the hydrophilic caps are constructed 
using the first bead of the DPPC head group (Q0, in the MARTINI topology). The 
protein beads also have bonded interactions where the bond length is 0.45 nm and 

the bond angle is set to 180°. The associated harmonic force constants for the bond 
_ 2 _ 2 _ 1 _ 2 

lengths and angles are 1250kJmol nm and 25kJmol rad . Based on the 
hydrophobic mismatch with the bilayers, the proteins are classified into three cate¬ 
gories: (I) positive mismatch (£ > T) Q ) (ii) negative mismatch (I < T) d) and (iii) 
no mismatch (£ ~ Ho)- To create different mismatches, we alter the number of 
beads in the protein core. These idealized proteins do not contain charges. 

Proteins are embedded in the equilibrated bilayer at 279 K. The resulting system 
is then heated to the required temperature and equilibrated for another 1.2/is. All 
the subsequent averages are performed using 10 independent trajectories each 600 ns 
long. The assembly of proteins is also performed using the same DPPC bilayer system 
with 3200 lipids and 50000 water beads. In this case, two proteins are inserted in this 
bilayer with centers at a distance of 14 nm and the simulation is carried out at 309 K. 

The flat interface is stabilized by juxtaposing an ordered bilayer equilibrated at 
285 K and zero lateral pressure with a disordered bilayer equilibrated at the same 
conditions corresponding to the cooling and heating curves of the hysteresis loop in 
Fig. SI, respectively. The system thus constructed is equilibrated in the ensemble 
with fixed temperature, volume and numbers of particles. This ensemble allows for 
maintaining an area per lipid intermediate between the two phases, thus stabilizing 
the interface. 

Instantaneous interface. For the purpose of obtaining a smooth and continu¬ 
ous interface, (/>(r) is coarse grained by replacing replacing Dirac’s delta function 
with a finite-width Gaussian, (1/27 t£ 2 ) exp (— |r| 2 /2£ 2 ) . The replacement 

changes (/>(r) to (/>(r). The coarse-graining width, £, is chosen to be the av¬ 
erage separation between tail-end particles l and j when ((( f>i — ~ 

~ (4>i)) 2 ) in the ordered phase is 1/10. This choice yields a value 
of £=1.5 nm. The instantaneous order-disorder interface is the set of points S satis¬ 
fying £)(s, t) — (0d + (f) o)/2. Here, (ftd and 4>o are (0(r)) evaluated in the 
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disordered and ordered phases, respectively. At zero lateral pressure and 294 K, we 
find (j )d = 0.4zb0.02nm 2 and (f) 0 = 2.15=b0.2nm 2 . For numerics, a square 
lattice tiles the average plane of the bilayer, and the coarse-grained field (/>(r) is 
evaluated at each lattice node. Values between are determined by interpolation. For 
convenience, the Gaussian function is truncated and shifted to zero at 3^. Any value 
of £ within the range, lnm < £ < 2nm gives nearly identical 0(r). Outside that 
range, larger values obscure detail by excessive smoothing, and smaller values obscure 
detail by capturing a high density of short-lived bubbles of disorder. 
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In this Supporting Information, we examine hysteresis in the bilayer system, show results for mean-square displacements 
as functions of time, and provide evidence that a small transmembrane ct-helix protein can be orderphobic if the side chains of 
the protein disrupt the packing of lipid tails. 


The lipid bilayer system exhibits hysteresis 

Fig. SI A shows the change in area per lipid with temperature while heating and cooling a bilayer. There are finite jumps in 
area per lipid as the system transitions between the two phases, suggesting a first-order phase transition. Hysteresis occurs 
because ordering from the metastable disordered phase is much slower than disordering from the metastable ordered phase. 
Due to the difference in time scales, when contrasting melting and freezing from heating and cooling runs, the melting points 
from heating runs as shown in Fig. 1A provide the more accurate estimates of the actual phase boundaries. Systematic errors 
due to small system size and heating rate have not been estimated. 
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Fig. SI. Structural measures of different phases as a function of temperature, T. (A) Variation in area per lipid with temperature during heating and cooling shows finite 
jumps and hysteresis. (B) Average local orientational order, (4>i), also shows finite jumps as a function of temperature while heating and cooling. Magnitudes of heating and 
cooling rates are 3 K/fis. 


The pre-melting layer has a higher mobility than the ordered phase 

In Fig. S2 we show the mean squared displacements of lipids in the bulk ordered phase, bulk disordered phase, and the 
pre-melting layer induced by an orderphobic protein. As discussed in the main text, an orderphobic protein increases the 
mobility of the lipids in its vicinity. Note that the center of mass of the membrane fluctuates in time. These fluctuations affect 
the absolute positions of lipid molecules, but they are irrelevant to the issue of lipid mobility. Therefore, the mean-square 
displacements considered in Fig. S2 are for tail-end particle positions relative to the instantaneous position of the membrane’s 
center of mass. That is to say, for (|E(t) — E(0)| 2 ), where r i(t) is the position at time t of the Ith tail-end particle less that of 
the membrane’s center of mass. 

For large enough times, t, the mean-square displacements are asymptotic to 4 Dt, where D is the self-diffusion constant. For 
the disordered liquid phase, we see from Fig. S2 that the asymptotic region is reached within 10 2 ns, and that D 4x 10 -7 cm 2 /s. 
In contrast, for the ordered phase, the diffusive asymptotic limit is not reached until 10 4 ns, and D « 2 x 10 _9 cm 2 /s. The 
mobility of lipids within the disordered layer surrounding the orderphobic protein is an order of magnitude larger than that of 
lipids beyond that region and in the ordered phase. 

Experimental results for lipid diffusion constants in disordered and ordered bilayers are D « 3 x 10 _8 cm 2 /s and D ^ 
2 x 10 _lo cm 2 /s, respectively [61]. The simulation is in harmony with experiment for the two-order of magnitude difference 
between the ordered- and disordered-phase values of D. Of course, absolute values of D are beyond the scope of what can be 
predicted from our simulations because coarse graining omits degrees of freedom that would increase friction and decrease D. 

A model Cf-helix is orderphobic 

Here, we consider the MARTINI model for KALP23 — a polypeptide chain with 23 residues, consisting of alanines and leucines 
flanked by lysines [62]. This molecule has a hydrophobic thickness of i ~ 3.0 nm, which means that it has essentially no 
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Fig. S2. Mean-square displacements as functions of time, t, for lipids in the disordered phase, the protein induced disordered domain, and the ordered phase. The functions 
are shown on log-log scale (main graphs) and linear scales (inset). 


hydrophobic-length mismatch with the ordered bilayer. Nevertheless, it is orderphobic because its side chains perturb the 
ordered lipid phase to an extent that a pre-melting layer is formed around the protein. This behavior is demonstrated with 
the aid of Fig. S3. The panels render configurations from a simulation in which we have placed this model protein in the 
MARTINI model for the ordered DPPC bilayer system considered in the main text. The pre-melting layer that forms around 
the protein causes the protein to tilt so as to keep its full hydrophobic length in contact with hydrophobic tails of the lipids. 
The interface separating its disordered domain from the surrounding ordered phase remains stable throughout a molecular 
dynamics trajectory running for more than 1 /xs. 



Fig. S3. A model o-helix is orderphobic. (A) Cross section of ordered phase of a hydrated DPPC membrane containing one transmembrane KALP23 protein. Solvent water 
is not rendered for purpose of clarity. Configuration was obtained after running simulation for roughly 1 fis at 294 K and zero lateral pressure. (B) Configuration of tail-end 
particles for the top monolayer, with grey points locating the instantaneous interface. 
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